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Abstract 

We propose classical equations of motion for a charged particle with 
magnetic moment, taking radiation reaction into account. This generalizes 
the Landau-Lifshitz equations for the spinless case. In the special case of 
spin-polarized motion in a constant magnetic field (synchrotron motion) we 
verify that the particle does lose energy. Previous proposals did not predict 
dissipation of energy and also suffered from runaway solutions analogous to 
those of the Lorentz-Dirac equations of motion. 



1 Introduction 

What is the classical motion of an electron in a constant magnetic field? Surpris- 
ingly, this question does not have a definitive theoretical answer even now. The 
complications arise when the radiation from the particle is considered. Without 
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it, the answer is quite elementary: a helical path around the magnetic field[Q]| 
obtained by solving the Lorentz equation of motion^ 



du^ q v 
dx m 



-Fuvu v (1) 



But we know that matters are not so simple in reality: an accelerated charge radi- 
ates, and the radiation emitted by it exerts a force back on it. If the energy is not 
replenished the path should be a spiral, the radius of the circle decreasing with the 
energy. 

Determining the force exerted on a point particle by its own field is fraught 
with conceptual difficulties: the force is infinite. Dirac[2] made a major step 
forward by showing that this divergence can be removed by a renormalization of 
the mass. The resulting equation (the Lorentz-Dirac equation of mo tion)| 



duy. _q , } d 2 u y 2q 2 



is still not the correct answer. It predicts runaway solutions: even in the absence 
of any external fields, a particle can accelerate away to infinity with exponentially 
growing energy. Part of the problem is that the radiative terms persist even in 
the absence of external forces. Also, the equation involves the third derivative of 
position; so it needs a new initial condition. Dirac suggested that this be chosen 
such that the energy remains finite forever. However, such solutions suffer from 
another problem: a violation of causality; the particle would accelerate before the 
external field is turned on. 

The textbook of Landau and Lifshitz[3J proposed (without proof) a way out of 
this dilemna. Replace the derivative of acceleration by its leading approximation 
in an expansion in powers of e: 



m p 



2q 2 

e = ~- (3) 
3 m 



dx m 1 J dx 

The radiative term now vanishes when the extrenal force is zero; and the equa- 
tion now involves only the second order derivative of position. The free particle 
does not run away. Further analysis shows that the orbits in a constant magnetic 



'in units with c = 1. q and m are the charge and mass respectively. Also, u} 1 = is the four- 
velocity and t is proper time. Note the constraint 7] /iv M fl w v = 1 and its derivative *7/iv M ' 1£ yf = 0. 

2 The last term includes a projection operator that ensures that the constraint TJ/jv"' 1 ^- = is 
satisfied. 
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field [4] and in a Coulomb field [5] decay as expected. But why should we stop at 
first order in the iteration in powers of e? Although it works, the ad hoc prescrip- 
tion of Landau and Lifshitz is unsatisfactory. 

To get sensible solutions to the Lorentz-Dirac equation, its new degree of free- 
dom must be continuously fine-tuned so that the trajectory stays on a submanifold 
of the phase space (the critical submanifold). Spohn used modern ideas in the 
theory of dynamical systems to determine the equation of motion projected to this 
critical submanifold: it is precisely the above Landau-Lifshitz equation! Thus 
we consider the century-old problem of determining the equation of motion of a 
radiating point charge to be settled finally. 

However, the electron has a magnetic moment as well as electric charge; ra- 
diation arises not only from the acceleration of the charge but also the precession 
of the magnetic moment. In the relativistic case, these are of comparable magni- 
tude. Therefore, it is necessary to extend the Lorentz-Dirac and Landau-Lifshitz 
equations to the case of a spinning charged particle. 

Bargman, Michel and Telegdi[6] proposed, more than a half century ago, an 
equation for a (non-radiating) classical particle with charge and magnetic moment. 
However, it does not seem to be the classical limit of the Dirac wave equation of 
the electron. The correct equation must have a canonical form whose quantization 
leads to the Dirac wave equation. The BMT equation does not appear to satisfy 
this condition. 

A better approach was proposed more recently by Barut and Zanghi[7]. We 
will review their ideas, which show that a canonical quantization does lead to the 
Dirac equation. We will also find solutions of the Barut-Zanghi equation for a 
non-radiating charged particle in a constant magnetic field: itself a new result. 

Barut and Unal[[8l also derived the analogue of the Lorentz-Dirac equation, 
after a mass renormalization. However, it suffers from the same runaway solutions 
of the free particle, as the spinless case. 

In this paper, we propose an equation of motion for a radiating spinning charged 
particle with gyromagnetic ratio 2 (e.g., the electron or the muon) by applying a 
prescription analogous to that of Landau and Lifshitz to the Barut-Zanghi equa- 
tion. We will show that, at least for spin-polarized orbits in a constant magnetic 
field, the orbits decay as expected physically. It should be possible to test our 
equations using a small (few MeV) electron synchrotron or a muon storage ring, 
such as the one used for the determination of g — 2. Also, it should not be difficult 
to revise the predictions for an electron accelerated by a laser fields. Current 
predictions do not include the effect of the magnetic moment. 

We do not attempt, in this paper, to derive our prescription from first principles 
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as Spohn did. 



2 Non-Radiating Spinning Charged Particle 

Can we talk of the motion of a particle with spin in the classical limit? At first this 
would appear to be impossible. But we do it all the time. The motion of electrons 
in synchrotron of a radius of several meters (if not kilometers) is surely classical: 
even if the spin were polarized, the differences in the quantum energy levels of 
the electron are much too small to be of significance. If we can find a lagrangian 
whose canonical quantization is the Dirac wave equation, it can be the basis for 
such a classical treatment. 

The lagrangian proposed by Barut and ZanghHis, 

L=^[zz-&]+Pp^-H, H = zfz\p v -qA v ], z = Z f % (4) 

Here z is a Dirac spinor with four complex (not Grassmann) numbers as com- 
ponents. The dimensions of z are that of the square root of action (=mass* length) 
and that of x is the inverse of mass, in units where c = 1 but fi ^ 1 . Then qF^ v has 
dimensions of mass over length. 

The equations of motion are, with itn = — qA^ 



zj^z 



dx 

.dz n 

-si = f "" z 

= qF, v ZYZ 

(5) 

Note that, x is not proper time as is not necessarily of constant length. But 
we have the orthogonality condition 

x% = 0. (6) 
These equations follow from the Hamiltonian 

H = zfz%n (7) 



They use a slightly different system of units. 
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and Poisson Brackets (P.B.) are obtained using 



{n^,n v } = -qF^ V: {n^,x v } = 8^ 

{z,z} = i, {^,z}=0 = {z,x} (9) 
All other pairs of P.B. are zero. We can introduce a constraint 

zz = a (10) 

consistent with the P.B. and the equations of motion. The quantity a, with 
the dimensions of action, sets the scale for the intrinsic angular momentum of the 
particle. 

Upon quantization, the wave function can be thought of as a function of x, 
Pli = it is also a polynomial in z, with z = — hjj= satisfying the constraint 

d 

-hz^-y = aw. (11) 
oz 

Thus the parameter a is quantized in multiples of h. The smallest value is 
a = h for which yf(x,z) = zw(x) is a linear function. 

The Dirac wave equation is the eigenvalue equation for H, with eigenvalue 

am: 

ihy 11 - — w = m\if (12) 
dxV 

Other quantizations, for which yr(z) is a higher degree polynomial, correspond 
to the wave equations of Bhabha and Harish-Chandra. They could be interesting 
as wave equations for composite particles like hadrons[10]. 

It is convenient to introduce the velocity and spin variables 

v »=zfz, S» v = l -z[f,f]z (13) 

so that 
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AS^Tty 

qF j _ lv v v 

(14) 

These equations describe a relativistic analogue of the isotropic precessing top 
of Poinsot ifTTTl . as we explain in the Appendix. The hamiltonian 

# = v% (15) 

describes a coupling to the external field which causes a precession of the 
velocity and spin. For a constant electromagnetic field, K = % ■ % — qS^ v F^ v is 
conserved as well: 

^- [n-n-qS^F^] =2qF^v v -2qF^ v v v n^ = 0. (16) 

CI L 

This is the classical analogue of the square of the Dirac operator, correspond- 
ing to the gyromagnetic ratio 2. We will see later that it is proportional to the 
angular momentum in the plane of the electromagnetic field. 



dx 

dx 
diz^ 

dx 
dS» v 

dx 



3 A Non-Radiating Spinning Particle in a Constant 
Magnetic Field 

Now let us turn to solving the equations of motion in the case of a magnetic 
field pointed along the third direction: F\2 = —F21 > 0; and all other components 
are zero. Before we solve the problem including radiation damping, we must 
recover the solution without radiation. Even this appears to be new: Barut and 
collaborators only solved the case of a free particle. 

We consider only orbits that lie in the 12 plane. That is, we restrict to the 
case where all the tensor indices /x, v,p • • • are in the 012 subspace. It is useful to 
introduce the 'dual' variables 
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5M v 



— £[ivpF p 
= e^ vp S n 



Sn = 



(17) 



For a constant magnetic field, only the time component F° is non-zero. We 
can assume it to be positive without loss of generality. 
Then 



If we define 



they reduce to 



dit. 



p 



dx 

dv^ 

~~dx 
dx 



£^v P v v FP 



dlty, 

dx 
dX» 

dx 
dpp_ 
dx 



-2£» vp X v n p 



2e }1V pp v K p 



(18) 



(19) 



3.1 Spin-Polarized Orbits 

There is a special class of 'left-handed' orbits for which = 0. (Also, the op- 
posite case with = 0.) For these, the spin and velocity are synchronized in the 
following way: 
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V = \^v P v p (20) 

In the non-relativistic limit, this means the spin is pointed along the magnetic 
field (times q). 

Since these solutions are periodic (not chaotic), orbits of this kind in a syn- 
chrotron should be of special interest. If the gyromagnetic ratio is slightly differ- 
ent from two, the spin would depart from this orientation, which is the principle 
behind a famous experiment for the measurement of g — 2 for the muon. 

The equations reduce to 



dK » _ P v v F p 

^ = -2e^ v Pv v K p 
dl y 

(21) 

The tensor e^ v = e^vprfP 17 form the structure constants of the three dimen- 
sional Lorentz Lie algebra 50(1,2). This is analogous to the way that the Levi- 
Civita tensor gives the stucture constants of the Lie algebra of rotations in Eu- 
clidean three dimensional space. Denote the Lie product of SO( 1 , 2) as a relativis- 
tic cross product 

(w x w) M = £^ p u v wp (22) 

This differs from the familiar cross product in Euclidean geometry by some 
signs: 

{uXw)q = U\W2 — U2W\, (uxw) l = —U2W0 + UQW2, (u X w) 2 = ~ UqWi + U\ Wq 

(23) 

Then our equations are a relativistic version of the equation for a precessing 

top: 

t = v x F, v = -2v x % (24) 

It follows that 

vv = (25) 
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so that v • v = a 2 is conserved; a is some number (like Planck's constant) with 
the dimensions of action; the spin is of magnitude |. Also, the parameter x is 
proportional to proper time in this special case. 

The energy 7Zq = E and 

K=K il Tl^-v^F ll (26) 

are conserved^. Now, 

vo = 2(v 2 ^i -v\K 2 ) (27) 

and 

vI = a([v\ + v 1 1 \ + -[vi^+v 2 ^ 2 ] 2 ) (28) 
Re-expressing the r.h.s. in terms of conserved quantities, 

vo = 2 yf [v 2 - a 2 ] [E 2 -K + v F ] -[H- Ev Q f (29) 

Weierstrass elliptic function satisfy this differential equation. The remaining 
components of n,v can be obtained by solving the other constraint equations. It 
is also possible to understand the solution as an integrable hamiltonian system, a 
relativistic analogue of the precessing top. (See the Appendix) 



4 The Equation of Motion for Radiating Spinning 
Particles 

Barut and Unal followed the method of Dirac to derive the analogue of the Lorentz- 
Dirac equation of motion including the self-force. Only the equation for %» is 
affected: 

dftii v ? r v u v v^ | 2 v v 9 v-vv v I 

-7T=^mW +<? 2 T7MV-— \~ 7/ <2( (30) 

ax i v-vJ I 3v-v 4(yv) J 

The last term is the contribution from self-interactions, after a renormaliza- 
tion. It is a singular perturbation i.e. q 2 is small and it changes the order of the 
equations. It has runaway solutions, just like the Lorentz-Dirac equation. Let us 

4 K is proportional to angular momentum. 
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apply to it a prescription like that of Landau and Lifshitz, to get a second order 
equation of motion. 

Naively, the Landau-Lifshitz prescription amounts to using 

v^=4S^n v + 4S^7t v (31) 
and replacing % by its zeroth order contribution: 



v'M 4S^ v 7t v + AS^qFypvP = 4 [v v 7^ - v^^} K v + AS^ v qF vp v p (32) 

But this can't be right: for a free particle {F^ v = 0), the radiation reaction force 
should be zero. Dropping the terms not proportional to F^ v , 

^ = q F, vV v + k v - ^1 SW pa v° (33) 
ax 3v- v L v- v J 

Along with the previous equations for v and 5, this system of ordinary differ- 
ential equations is our proposal for the equation of motion of a spinning radiating 
charged particle. We must verify that the solutions are physically sensible: there 
should not be runaway solutions. We will study the case of polarized orbits in a 
constant magnetic field and show that the energy decreases with time. The ulti- 
mate test of our proposal must be experimental, especially since we do not yet 
have a derivation from first principles. 



5 Radiating Spin-Polarized Particles in a Magnetic 
Field 

Now let us turn to our proposed equations for a radiating spinning charged parti- 
cle. Again, we will look for solutions in the special case where the spin is polar- 
ized 



S» = (34) 



1 

2' 

Since the radiative terms do not change the equations for , S^ v this condition 
is still preserved by time evolution i.e. 5^ = ^v^. Unpolarized orbits are chaotic 
and the radiation damping is likely to be even larger. (See Appendix) 

The equations then reduce to 
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ax 3v ■ v L v • v J 



We know that, v ■ v = and v ■ v = Vq — (v\ + v 2 ,) = a 2 



(35) 



^ = e M v P v v F p + [a 2 T7 MV - v M v v ] yF v (36) 



where 



_ 4q 2 
Y ~3a 2 

In terms of relativistic cross products 



(37) 



fc = vx [F + yFx v], v = -2v x tz; (38) 

Thus, even the radiation damping terms have a simple meaning in terms of the 
relativistic cross product of the Lie algebra SO (1,2). 
Note that 

H = vk (39) 

is still conserved and m = ^ continues to have the physical meaning of mass. 
But energy E = Kq and K are no longer conserved. The evolution of vo continues 
to be determined by the equation we derived earlier. Thus we get the system 



dE , 9 9n 

- = —yFo(vQ — a) 

^ = 2yF [a 2 E-Hv ] 

v = 2sJ[v 2 - a 2 ] [E 2 -K + v F ] -[H- Ev Q ] 2 

(40) 

Since E and K are no longer constants, the solution is no longer an elliptic 
function. Also, noting that Fq > in our convention and v\ — a 2 = (v\ + v\) > 0, 
we see that energy is monotonically decreasing. 
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Appendix: Hamiltonian Formalism 

The Barut-Zanghi equations for a spinning charges particle follow from the P. B. 
forming the Lie algebra 50(2, 3) : 



{vV v } = 45^ v 
{S^\S pa } = n^S va + ri va S^-ri v PS^ a -'n^ a S v P 

(41) 

Using previous relations we can also find the following P.B. : 



{v M ,v v } = 45 MV 
{5^,5y} = 
{S M ,v v } = e^ivav a 

(42) 

which generalizes the rotation algebra. When the third axis is ignored, this 
reduces to 50(2,2) 50(1,2) ©50(1,2). This splitting can be seen by in terms 
of the linear combinations , p M : 



We have also the P.B., 



{A^, Ay} 

{Pm>Pv} 

{A M ,pv} 



^/xvctA 17 




(43) 
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{A M ,7Tv} = 
{p M ,%} = 

{n^itv} = -£^ vp F p 

(44) 

The hamiltonian is then 

H = 2[X-p]-n (45) 
Non-Radiating Spin-Polarized Orbits in Constant Magnetic Field 

Since 

v° = y/a 2 + v b v b (46) 

and Kq = E is conserved we only have four co-ordinates in the phase space: 
7t a ,v a . The P.B. are 

{K h K 2 } = -F (47) 

and 

{vi,v 2 } = 2y / a2 + v 2 +v 2 (48) 

The latter corresponds to the standard symplectic form on hyperbolic space. It 
is not hard to find canonical co-ordinates 

{-7=) -7=} ~~ ^ { a/2/? cos 0, y / 2psin0},{vi,V2} — > {a sinhO cos0,a sinhO sin^} 



so that- 

{p,0}=l, P = ^ F (rf + n%), = arctan^ (49) 

2Fq Tti 

{vo,(j)} = 1, = arctan — (50) 

vi 

Note that J = p + ^vq = ^ is the generator of rotations; it is a conserved 
quantity corresponding to angular momentum. The hamiltonian is then 
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H = Ev Q -\j 2F [v 2 Q - a 2 ] (7 - - v ) cos ( 9 - 0) 

Since 

pd6 + lp?od<l) = Jd6 +voda, a = ^(0-0) 
we can use the pairs 7, and vo, oc as canonical variables: 

H = Ev - yj 2F [v§ - a 2 ] (7 - ^v ) cos 2a 
The equation of motion is 



Thus 



= 4F [v§ - a 2 ] (27 - Vo ) - 4 (tf - Ev y 



(51) 



(52) 



(53) 



(54) 



(55) 



which is the Weierstrass equation. Because of its polarization, this also deter- 
mines the spin. 
Also 



Thus 



dd _ dH 
dz~~dJ 



iFo 



v o~ a 



27 — vq 



cos 2a 



dG H-Evq 



(56) 



(57) 



dx (27- vo) 

which determines also in terms of elliptic integrals. 
To determine the position variable, we note the that 

^[^-£nv P x v FP] = (58) 
This leads to the conserved quantity (corresponding to translation invariance) 

X il = K ll -£ llvp x v FP (59) 
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Thus (in a co-ordinate system in which Fq > 0,F a = and X c = 0) we have 



x b = -^r£ bc 7Z c (60) 

In other words, knowing ii\ , % is the same knowing the spatial orbit, except 
for a rotation by 90° and a scaling by Fq. In polar co-ordinates this amounts to 
knowing p = J — jVq and as functions of T, a kind of elliptic curve: 



= 2F [4(J-p) 2 -a 2 ]p-(H-2EJ + 2Ep) 2 , d ± = *LJL + | 

(61) 

It is not correct to interpret the x as the position of the electron: due to the 
phenomenon of zitterbewegung. Even for a free particle, x^{x) is a spiral and not 
a straightline. Barut et. al. have shown how to recover a center of mass variable 
for the electron from this complicated solution. 



A Symplectic Reduction of the General Planar Orbit 

The general case of planar but not spin-polarized orbits does not appear to be 
integrable. It corresponds to chaotic motion of an electron in a synchrotron. We 
can use the conserved quantities 

tzq = E, X ■ X = a 2 , p -p = b 2 (62) 

to reduce it to a conservative system with three degrees of freedom: the phase 
space is a plane (n a ) times a pair of hyperboloids parametrized by (A^Pjx). It is 
easy to check that 



{ P ,e} = i (63) 

with 

1 7t2 
p = —K a K a , = arctan — (64) 
2F K\ 

Also, 
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{Ao,0} = l, <j> = arctan — 

M 
P2 

{pO,X} = h X = arctan — 

Pi 



(65) 



In terms of these 



H = 2E(Xo-po)-2(Xa-pa)n a 



= E(Xq - p ) - yj2F Q [X 2 - a 2 ] pcos[0 - 0} + ^2F [p 2 - b 2 } pcos[ X ~ 0} 

(66) 

Another conserved quantity is the generator of rotations^ 

J = p + Xo + p (67) 

The reduced system has three degrees of freedom and two conserved quantities 
H,J : just one conserved quantity short of an integrable system. When Xq = a 
(or po = b ) one degree of freedom decouples and we get an integrable system. 
Circular orbits correspond to the case = %, Xq + po = constant. 

It may be more convenient to use a = — , /3 = % — 6,6 as the independent 
position variables. Then 

Pada + ppdfl+pedd = pdd + X^dfy + p dx = pdO + Xod[0 + a] + pod[0 + j8] 

(68) 

Pa = Xq 
Pp = Po 

pe = p + Xo + po = J 

(69) 



5 This is linearly related to the quantity K we introduced earlier: 

K = 2FqJ 
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Then we get a system with two degrees of freedom with hamiltonian 

H = E{Xq- Pq )-^2F q [X 2 - a 2 ] [J - Ao - p ] cos a + J 2F [p 2 - b 2 ] [J-Xq- p ] cos j8 

(70) 

This appears to be a chaotic system. In general, the information about the 
polarization of spin in the initial condition will be lost rapidly and the particle 
beam will look unpolarized. For the special case po = = b the spin remains 
polarized and solutions can be obtained in terms of elliptic functions as we saw 
above. 
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